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INTRODUCTION 


This  report  focuses  on  five  individual  work  areas:  (i)  synthesis  of 
S10x,  (1i)  synthesis  of  SiC,  (iii)  synthesis  of  Pd2Si,  (iv)  electrical 
diagnostics  of  films  and  (v)  computer  modeling  of  film  synthesis. 

Reactive  sputter  deposition  of  refractory  oxides,  under  conditions  in 
which  the  reaction  between  the  sputtered  species  and  the  oxygen  is  isolated 
to  the  deposition  surface,  has  been  experimentally  Investigated.  Results 
for  the  formation  of  S10x(0  <  x  <  2)  under  such  experimental  conditions 
are  discussed  below.  It  Is  shown  that  the  dependence  of  the  deposited  film 
stoichiometry  on  oxygen  pressure  and  the  incidence  rate  of  silicon  can  be 
correlated  with  reported  results  of  oxygen  adsorption  onto  single  crystal 
silicon  surfaces  and  amorphous  silicon  surfaces.  This  correlation  is  made 
In  terms  of  oxygen  exposure,  t(Torr-sec).  For  our  work,  a  dynamic 
exposure,  Lq,  is  defined  for  steady  state  synthesis  conditions  by  placing 
the  origin  of  a  moving  reference  frame  at  the  deposition  surface.  By  means 
of  this  correlation.  It  Is  found  that  the  sticking  coefficient  of  oxygen 
onto  room  temperature  substrates  Is  dependent  of  oxygen  coverage  (i.e., 
surface  stoichiometry)  for  coverages  less  than  0.8.  This  requires  that  a 
mobile  precursor  adsorption  process  be  operative  for  oxygen  In  our 
experiments. 

The  reactive  sputtering  study  was  briefly  extended  to  an  N2 
environment  and  conditions  for  the  formation  of  SI3N4  determined.  In  the 
SIC  work,  the  emphasis  was  on  high  substrate  temperature  depositions  where 
surface  contamination  and  surface  conversion  of  Si  to  SIC  was  observed. 
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Much  effort  was  spent  on  high  temperature  development  and  on  maintenance  of 
a  clean  environment  in  this  high  temperature  domain.  Steady,  rather  than 
dramatic,  progress  is  being  made. 

Using  the  phase  spread  technique  and  a  room  temperature  substrate,  the 
structural  character  of  the  films  as  a  function  of  composition  in  the  Pd- 
Si  system  have  been  determined.  In  the  electrical  diagnostics  of  films 
area,  an  analysis  has  been  made  of  the  proposed  measurement  technique  for 
the  specific  case  of  in-situ  film  monitoring  during  oxidation.  The  results 
indicate  that  the  technique  Is  only  marginal  with  respect  to  sensitivity 
revealing  interface  capacitance  and  interface  resistance  effects. 

In  the  theoretical  section,  an  important  new  formulation  of  the  stress 
and  elastic  constant  tensors  in  the  vicinity  of  interface  is  given.  This 
is  applied  to  the  simple  example  of  an  absorbed  layer/substrate  interaction 
system.  For  the  Si-0  system,  radial  distribution  function  and  angular 
distribution  data  have  been  generated  for  bulk  a-quartz  and 
a-crlstoballte  while  surface  reconstruction  producing  peroxide  bond 
formation  has  been  found  for  the  (100)  a-cristobalite  surface. 


A.  Experimental  Results 

1.  Reactive  Sputter  Deposition  of  Oxide  and  Nitride  Films 
SI 0X  -  The  Silicon-Oxygen  System  Analysis 

Reactive  sputter  deposition  of  refractory  oxides,  under  conditions 
In  which  the  reaction  between  the  sputtered  species  and  the  oxygen  is 
Isolated  to  the  deposition  surface,  has  been  experimentally  demonstrated 
and  reported  on  In  many  previous  reports^.  The  experimental  results 
presented  demonstrate  that  the  synthesis  of  silicon-oxygen  films  under 
controlled  conditions,  designed  so  that  reaction  between  silicon  and  oxygen 


occurs  at  the  deposition  surface,  allows  formation  of  films  of  reproducible 
stoichiometry  and  high  uniformity  onto  room  temperature  substrates.  To 
date,  we  have  discussed  these  results  in  terms  of  the  relative  rates  of 
Incidence  of  the  two  reacting  species  onto  the  deposition  surface.  It  is 
Important  at  this  point  to  consider  experimental  results  reported  by  other 
Investigators  concerning  the  adsorption  of  oxygen  onto  surfaces  of  bulk 

single  crystals  and  of  amorphous  silicon.  Such  surfaces  are  static  in  that 
the  number  of  Si  atoms  is  fixed  by  the  character  of  the  initial  surface. 
In  our  experiments,  a  dynamic  moving  surface  is  studied  in  which  both 
silicon  and  oxygen  are  supplied  to  a  moving  surface,  i.e.,  a  film  is  being 
grown. 

More  sophisticated  experiments  in  which  ultra  high  vacuum  studies  of 
oxygen  absorption  on  high  index  silicon  planes  (2,3,4)  and  amorphous 

silicon^5)  have  been  reported.  In  these  studies  the  oxygen  is  observed  to 
be  absorbed  with  the  coverage  being  a  linear  function  of  oxygen  exposure  in 
Langmuirs  (torr-sec).  Sticking  probabilities  are  small  (2  x  10*4  to  3  x 
10*3  and  are  strongly  affected  by  the  presence  of  excited  oxygen. 

Additionally,  the  oxygen  is  absorbed  in  an  unreacted  state. 

This  reaction  (Si  +  02  ♦  Si 02 )  is  apparently  strongly  catalyzed^4)  by 
the  presence  of  small  amounts  of  carbon.  Adsorption  on  amorphous  silicon 
does  not  appear  to  be  substantially  different  from  absorption  on  (111)  or 

(110)  silicon  surfaces.  The  sticking  coefficient  of  oxygen  on  silicon 

(111)  and  (110)  surfaces  were  reported  to  be  2  x  10’4  and  6  x  10"3 
respectively.  The  higher  value  on  (110)  surface  can  be  explained  by  Its 
higher  Index  and  the  possibility  of  surface  step  effects. 

In  this  work ( 12 )  the  sticking  coefficient,  S,  was  determined  by 

analysis  of  experimental  data  obtained  using  Auger  Electron  Spectroscopy^* 
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(AES),  High  Resolution  Electron^2*2)  Spectroscopy  (HRES),  and 
Elllpsometry^2*3)  (E).  It  was  assumed  in  all  these  experiments  that  the 
coverage,  &,  was  directly  proportional  to  the  intensity  of  the  signal 
monitored  which  was  determined  as  a  function  of  oxygen  exposure,  L,  and 
that  unit  coverage  was  produced  when  the  experimental  signal  saturated. 
Under  these  assumptions  the  coverage  dependence  of  the  sticking  coefficient 
S (6 )  could  be  derived  as 


d  vp 

—  (0)  =  S{6)  —  (1) 

dt  N 

o 

on 

where  8  is  the  fractional  coverage,  P  the  pressure,  v  =  3.48  *  10 
molecule/cm2  -  sec  Torr  and  N0  the  number  of  surface  sites  available  for 
absorption. 

The  relationship  between  coverage,  8,  and  exposure  is  similar  to  the 
relationship  between  composition  and  relative  rates  of  incidence  shown  in 
Fig.  1,  A  more  direct  comparison  can  be  made  if  our  data  is  displayed  in 
terms  of  the  oxygen  concentration  given  as  X  in  Si 0X  with  0  <  X  <  2  and  a 
dynamic  exposure  Lq.  The  correspondence  between  coverage,  6,  and  X  in 
S10x  Is  straight-forward  If  we  assume  that  the  absorbed  species  is  02  and 
that  unit  coverage  Is  attained  when  there  Is  one  oxygen  molecule  absorbed 
for  each  silicon  surface  atom. 

Dynamic  exposure,  Lq,  Is  defined  by  the  time  necessary  to  deposit  a 
monolayer  of  silicon  atoms  which  is  used  In  the  calculation  of  oxygen 
exposure  at  pressure  P.  The  number  of  atoms  on  an  exposed  (111)  silicon 
surface  Is  calculated  to  be  ~  3.96  *  1014/cm2.  The  rate  of  incidence  of 
silicon  atoms,  v$1,  has  been  already  given  as  v$.  =  R$. /v$i  where  RSi 
Is  the  deposition  rate  In  A^/sec  and  and  Is  atomic  volume  of 


silicon.  The  time  t,  to  deposit  one  monolayer,  i.e.  3.96  x  10*4  atoms/cm^ 


is  then 


3.96  *  1014  3.96  *  1014V$. 


Thus,  a  dynamic  exposure  Lq  for  oxygen  can  be  defined  as 


LD  =  po2  T 


where  Pn  is  given  in  torr  and  t  in  seconds.  During  a  steady  state 
u2 

synthesis  experiment  static  exposure  and  dynamic  exposure  are  analytically 
equivalent  at  small  time  steps. 

our  experimental  results  are  plotted  as  X  In  Si 0X  versus  dynamic 
exposure  Lq  in  Fig.  2.  The  results  of  Ritter^8)  are  also  plotted  for 
comparison.  The  differences  are  striking.  First,  the  dynamic  exposure 
range,  in  Langmuirs,  over  which  our  data  extends  is  approximately  10-4  to 
10-8  L.  In  order  to  gain  the  same  results,  Ritter's  range  was 
approximately  10~8L  to  10"PL.  No  clear  explanation  for  this  appears  in  the 
data  of  either  Investigation.  It  is  likely,  though,  that  Ritter's^6)  work 
had  substantial  amounts  of  excited  oxygen  present.  The  absorption  and 
reaction  probability  of  excited  Og  Is  orders  of  magnitude  higher  than 
unexcited  Og.  Therefore,  the  reactivity  of  the  oxygen  in  the  evaporation 
experiment  was  undoubtedly  strongly  Increased  by  the  experimental 
arrangement  used.  In  our  work,  the  oxygen  was  introduced  at  the  substrate 
In  an  unexcited  state.  The  substrate  to  sputter  source  spacing  was  large 
enough  that  the  coupling  between  the  plasma  and  the  deposition  surface  was 
minimal  (a  characteristic  of  planar  magnetron  systems).  Additionally, 
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substrate  heating,  consistent  with  the  deposition  parameters  (Ts<  65°C)  was 
observed  for  extended  synthesis  experiments  (<  2  hrs)  substantiating  that 
the  source-substrate  coupling  was  minimal. 

It  is  informative  to  compare  coverage^2,3)  as  determined  by  E,  HRES 
and  AES  with  our  data  for  X  in  Si 0X  versus  exposure,  LD,  as  shown  in  Fig. 
3.  The  qualitative  behaviors  are  very  similar  and  quantitative  agreement 
occurs  with  respect  to  the  exposure  range  involved  (i.e.,  10“4  to  10~3 
torr-sec).  This  implies  that  the  rate  controlling  step  in  the  synthesis 
process  is  identical  to  that  in  the  absorption  process  of  02  in  silicon. 
Hence,  our  rate  controlling  step  Is  the  absorption  of  02  with  reaction  to 
form  S10x  occurring  once  the  absorption  step  has  taken  place. 

The  slopes  of  the  film  stoichiometry  curves  (X  in  Si 0X)  are  steeper 
than  the  absorption  curves  Implying  a  somewhat  larger  sticking 
coefficient.  When  plotted  In  a  linear  manner.  Fig.  4,  the  relationship 
between  x  In  S10x  1$  observed  to  be  a  linear  function  of  exposure  without 
experimental  scatter.  This  states  that  the  sticking  coefficient  of  oxygen 
Is  Independent  of  coverage,  a  condition  which  occurs  when  there  is  a  mobile 
precursor  to  the  che*.'. sorption  event.  Therefore,  It  is  quite  likely  that. 
In  our  work,  a  mobile  percursor  of  02  adsorbed  onto  the  surface  exists  and 
facilitates  the  chemisorption  process. 

A  sticking  coefficient,  S,  for  02  onto  our  moving  deposition  surface 

can  also  be  estimated  for  these  data.  If  we  assume  unit  coverage,  when  x  « 

2  for  x  In  SiOx,  the  coverage  can  be  linearly  related  to  x  by  3  *  x/2. 

At  constant  pressure,  P,  the  time  rate  of  change  of  the  coverage,  e,  can 

be  related  {Eq.  (1))  to  the  sticking  coefficient,  S(®).  The  sticking  co¬ 
efficient  derived  on  this  basis  Is  In  the  range  1.56  x  10"3  to  4.17  x  10-3, 
These  number  »  a  actor  of  5  to  20  larger  than  the  value  reported  by 


Ibach  (S  =  2.2  x  10"^)  for  an  atomically  flat  (111)  surface.  This  increase 
in  S  may  be  due  either  to  the  presence  in  our  work  of  activated  oxygen  or 
to  an  enhanced  reactivity  of  our  amorphous  Si 0X  surface  relative  to  ideal 
crystalline  surfaces.  The  latter  of  these  explanations  is  consistent  with 
the  effects  of  surface  steps  on  the  sticking  probability.  It  would  require 
6-8%  of  the  surface  atoms  to  be  associated  with  surface  steps  to  achieve 
sticking  probabilities  as  we  have  Inferred.  Computer  simulation  studies  of 
amorphous  surfaces  by  Garofolini  et  al^7)  tend  to  support  this 
interpretation. 

The  observed  independence  of  the  sticking  coefficient  on  6,  for  02 

on  the  silicon  surface  can  be  explained  by  a  mechanism  in  which  a 
physisorbed  mobile  presursor  state  for  oxygen  exists  on  the  sample  surface 
prior  to  the  chemisorption.  This  allows  02  molecules  to  move  on  the 

surface  so  that,  at  substantial  coverages  (e  *  0.5),  the  remaining 

unoccupied  chemisorption  sites  are  accessible.  A  simple  model of  the 
combined  physi sorption-chemisorption  process  shows  the  sticking  coefficient 
S(0)  for  the  chemisorbed  stated  to  be  given  by 

u(l  -  6) 

S(e)  =  - ^  (4) 

(1  ♦  rd/rc)  -  9(1  +  rd/rm) 

where  0  is  the  coverage,  a  the  fraction  of  incident  gas  molecules 
entering  the  precursor  state,  rd  the  rate  of  desorption,  rm  the  rate  of 

migration,  and  rc  the  rate  of  transition  from  the  precursor  to  the 

chemisorbed  state.  In  the  limit  of  6  ♦  0,  the  sticking  coefficient  S0  is 
given  by 


S 


o 


a 

0  +  Vrc> 


(5) 
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The  relationship  between  and  rc  is  easily  explained  in  terms  of  a 
potential  energy  diagram  as  seen  in  Fig.  5.  In  this  figure,  the  potential 
energy  of  an  adsorbed  oxygen  molecule  is  plotted  as  a  function  of  the 
molecule's  distance  from  the  solid  surface.  In  Fig.  5,  qp  represents  the 
heat  of  physi sorption  and  Ep  the  activation  energy  for  transition  from  the 
physisorbed  to  the  chemisorbed  states.  If  we  assume  the  ratio  of  the  pre¬ 
exponential  factors  of  rd  and  rc  is  then 


rd/rc  =  Kie*p[-(qp-ep)/RT>] 


where  R  is  the  gas  constant  and  T  the  absolute  temperature.  Equation  (5) 
then  becomes 


0  (1  +  Kiexp(-(qp-Ep)/RT) 

In  all  likelihood  a  «  l  so  that  the  quantity  (l+K1exp(-(qp-Ep)/RT)  is 

_3 

equal  to  approximately  500  since  S  -  2  x  10  .  Therefore,  to  a  good 

approximation. 


Sn  =  -  exp[q  -E/RT] 


indicating  that  chemisorption  occurs  with  an  activation  energy  close  to 
zero  (l.e.,  qp  >  Ep),  a  phenomenon  found  in  the  oxidation  of  many 
metals^®"**). 

The  result  of  this  Interpretation  is  that  the  residence  time  of  a 
physisorbed  oxygen  molecule  on  the  surface  and  its  migration  rate  determine 
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the  number  of  chemisorption  sites  accessible  per  physisorbed  molecule  and 
therefore  the  chemisorption  rate.  Hence,  as  the  substrate  temperature  is 
increased,  at  fixed  oxygen  pressure  and  silicon  incidence  rate,  the 
stoichiometry  of  the  deposit  will  decrease  as  is  observed  in  the  reactive 
deposition  of  many  oxides  including  reactively  deposited  SiOx.  This  is  due 
to  the  decrease  in  residence  time  of  physisorbed  oxygen  mobility  with 
increasing  temperature. 

(b)  Oxide  and  Nitride  Film  Deposition 

The  new  gas  flow  controller/argon  ring  set  up  at  the  target 

location  was  tested  to  determine  Pn  (max)  before  target  oxidation 

u2 

occurs.  The  gas  flow  controller  was  found  to  be  stable  in  the  rf  fields 
present  but  some  small  indication  of  oxygen  contamination  in  the  Ar  line 
was  suggested  because  of  occasional  peaks  in  reflected  power  readings.  Si 
deposition  rate  tests  were  carried  out  both  with  and  without  (Ar  released 
at  the  chamber  wall  ~  6"  from  the  target)  the  Ar  ring.  The  without/with 
Ar  ring  data  were  4. 2/4. 3  A/sec  at  300  watts  (rf)  and  6. 2/5. 8  A/sec  at 
500  watts  (with  3  x  10-6  torr  base  vacuum,  4  um  maintained  Ar  pressure  and 
4"  source/substrate  distance).  The  reflected  power  was  generally  steady  in 
the  0-3.0  watt  range  except  at  500  W,  with  the  Ar  ring,  it  jumped  to  50 
W.  We  thus  concluded  that  02  was  occasionally  getting  into  the  Ar  line. 

In  order  to  work  easily  with  3"  diameter  wafers  and  to  extend  the 
substrate  temperature  range,  considerable  effort  was  expended  in  the 
development  of  a  new  substrate  heater.  Modifications  to  an  existing  two 
element  quartz  lamp  were  unsuccessful  but  the  development  of  a  carbon  cloth 
heater  was  fairly  successful.  However,  a  deposited  film  was  observed  on 
the  Si  substrates  placed  on  the  heater  during  bakeout  even  with  no 


sputtering  and  no  Intentional  gases  present.  Microprobe  detected  carbon 
and  silicon  in  the  film  while  EDAX  analysis  showed  no  obvious  metals  to  be 
present.  Based  upon  the  partial  success,  a  platinum  wire  heater  was 
designed  for  the  higher  temperature  work. 

Relatively  successful  attempts  have  recently  been  made  to  deposit 

Si 3n4  even  though  our  early  work  led  to  measured  values  of 

n  =  2.4  to  2.8  instead  of  2.03.  Using  250  W,  P^,r  »  4.0  urn  and 

PN  «  0.26  um  ,  source  to  substrate  distance  was  considered  as  the 
n2 

variable.  At  5",  no  index  could  be  measured  which  indicated  a  low  density 
Si-rich  film.  At  3",  the  deposition  rate  was  about  80  A/min  and  the  film 
index  (~  1200  A)  was  n  *  2.00.  At  4",  an  aluminum  wire  probe  was 
inserted  just  above  the  substrate  surface  to  provide  a  glow  discharge 
during  sputtering  as  an  aid  in  dissociating  the  N2.  In  this  case,  a 
deposition  rate  of  90  A/min  was  found  and  the  film  index  (~  1500  A) 
was  n  *  2.09.  Thus,  the  use  of  the  glow  discharge  was  important  since  the 
deposition  rate  increased  rather  than  decreased  by  9/16  in  going  from  3"  to 
4".  From  these  experiments  it  was  felt  that  nitride  films  from  silicon 
oxy-nitride  through  Si 3N4  to  Si-rich  nitrides  could  be  easily  deposited 
using  these  techniques  and  parameters. 

2.  Single  Source  Sputter  Deposition  of  SiC  Films 

For  high  temperature  operation  (T  >  900°C),  both  the  single  quartz 
lanp  heater  and  the  double  quartz  lamp  heater  were  found  to  be  inadequate 
as  a  substrate  heater  because  of  shorting  and  arcing  problems  during 
operation.  However,  after  considerable  development  effort  a  carbon-cloth 
heater  was  fabricated.  To  date,  it  has  proven  to  be  extremely  stable  in 
operation  up  to  900°C  in  a  plasma  environment.  The  present  temperature 


10 


limitation  is  a  matter  of  power  supply  capability  rather  than  heater 
capability. 


During  this  heater  development  study,  a  surface  contamination  problem 
was  noticed.  It  was  noticed  on  all  depositions  wherein  an  ion-gun  surface 
cleaned  substrate  was  used  above  T  =  600°C  and  on  substrates,  not  ion-gun 
cleaned,  heated  to  T  *  900°C.  In  the  latter  case,  a  substantial  surface 
conversion  of  the  Si  to  SiC  occurred  (~  2500  A  thick)  as  determined  by  x- 
rays  and  electron  microprobe  analysis.  This  SiC  layer  was  non-uniform  in 
thickness  and  triangularly  shaped  islands  ~  several  microns  in  size  were 
observed.  The  source  of  t  carbon  was  considered  to  be  organic 
contamination  in  the  chamber.  Less  contamination  was  observed  after  having 
scrubbed  the  chamber  and  the  masking  plus  heater  parts. 

A  Meissner  trap  was  made  from  coiled  stainless  steel  tubing  and 
mounted  on  the  Inside  of  the  chamber  near  the  outer  wall.  Flowing  liquid 
N2  through  the  coil  reduced  the  system  pressure  from  typically  2  x  10"6 
torr  to  5  x  10“7  torr.  Monitoring  with  the  gas  analyzer  (RGA)  showed  that 
the  trap  primarily  reduced  the  water  vapor  level  in  the  chamber.  This 
reduced  the  silicon  surface  conversion  problem  somewhat  but  not  completely. 


New  Depositions 

The  list  of  new  depositions  are  presented  in  Tables  I  and  II.  Surface 
conversion  of  Si  to  SiC  was  seen  to  an  even  greater  degree  at  900°C  than 
with  depositions  at  lower  temperatures.  A  methanol  rinse  seems  to  make  the 
depositions  less  crystalline,  perhaps  enhancing  the  surface  conversion. 
Sample  564,  with  an  HF  dip  surface  clean  and  a  900°C  SiC  deposition  without 
tin,  appeared  to  be  (via  x-ray  diffractometry)  as  high  quality  as  our  best 
earlier  deposition  (#82-233)  although  visually  the  film  was  not  uniform. 


No  surface  contamination  was  found  with  the  deposition  on  sapphire.  Via 
electron  microprobe,  tin  was  found  in  all  the  samples  597-601  (500°C  to 
800°C  depositions)  indicating  that  the  tin  target  power  of  ~  105  W  was 
too  high. 

From  Table  II  we  note  that  a  set  of  Si C  depositions  at  500  W  was  made 
on  a  variety  of  substrates.  The  thin  film  on  (0001)  sapphire  is  a  great 
Improvement  over  all  previous  depositions  according  to  both  the  high  x-ray 
count  and  a  Read  camera  pattern  which  showed  only  the  111  B  SiC 
reflection  in  a  tight  +  3°  distribution  about  the  substrate  surface 
normal.  Interestingly,  this  substrate  had  no  ion  gun  cleaning  prior  to 
deposition.  The  sheet  resistivity  value  for  82-124  was  found  to  be  quite 
low  at  2.7  *2  cm  compared  to  the  highest  p$  to  date  of  590  O  cm  (sample 
81-124  on  (lT02)  sapphire  at  700°C).  Sheet  resistivity  values  for  the 

depositions  on  oxidized  Si  were  found  to  represent  measurement  of  the 
substrate  rather  than  the  thin  SiC  film.  Finally,  the  deposition  on  spinel 
could  not  support  a  voltage  suggesting  a  serious  interface  problem.  This 
film  appears  to  be  highly  stressed  and  is  gradually  lifting  off  the 
substrate. 

A  number  of  SiC  depositions  that  were  examined  by  TEM  showed  a 
hexagonal  <*-SiC  component  not  seen  by  any  other  analysis  method  and  it  was 
thought  that  this  might  arise  from  the  extensive  ion  milling  involved  in 
sample  preparation  or  from  the  electron  beam.  A  number  of  very  thin, 
~  500  A,  depositions  of  SiC  on  Si  were  made  to  sort  this  out  (#619-623). 
No  TEM  analysis  has  yet  been  made  but  extensive  optical  microscopy  at  1550 
x  Indicated  distinctive  surface  features  in  the  micron  size  range.  The 
largest  features  were  triangular  and  appeared  brighter  in  the  microscope. 
Either  tin  or  a  different  phase  was  suspected  but  tin  was  not  detected  via 


dispersive  x-ray  analysis.  These  samples  look  interesting  and  may  provide 
Information  about  the  early  stages  of  growth  of  SiC  on  Si.  A  TEM  sample 
preparation  technique  without  using  ion  milling  is  underway. 


3.  Synthesis  of  Pd  Silicide  Films 

<1 1 1>  Si  wafers  were  scribed  and  inserted  in  the  sputtering  chamber  as 
described  in  the  previous  report  (*2).  The  wafer,  which  was  located  at  the 
point  where  Pd2Si  was  calculated  to  be  deposited,  was  etched  in  dilute  HF 
until  the  surface  was  hydrophilic  and  then  the  system  was  immediately 
sealed  for  pump  down.  Material  was  deposited  at  room  temperature  or  at 
400°C  for  93  minutes  to  produce  approximately  1  wm  of  Pd2Si  so  that 
electron  probe  microanalysis  could  eventually  be  carried  out  on  the  film. 

The  results  of  the  phase  spread  deposition,  as  determined  by 
microprobe  analysis,  agreed  within  experimental  error  with  the  calculated 
results.  The  compositions  shown  in  Table  III  were  analyzed  using  x-ray 
diffraction.  The  Pd  lattice  was  found  to  expand  linearly  at  ~  0.20%  per 

atomic  percent  of  added  Si  In  the  1  to  10  At  %  Si  range.  At  about  10%  Si, 

amorphous  material  appears.  This  is  expected  because  of  the  deep  eutectic 
at  ~  15  At  %  Si.  At  ~20  At  %  Si ,  Pd2S1  is  first  observed  and  this  phase 
Increases  in  Intensity  until,  at  27  At  %  Si,  the  material  is  entirely 
Pd2S1.  The  material  then  transforms  into  PdSI  at  ~49  At  %  Si  and  it  is 
unclear  as  to  the  exact  progression  of  phases  between  Pd2Si  and  PdSi  as  the 
Si  increases. 

Initially,  when  the  samples  had  been  freshly  prepared,  there  were  a 
great  number  of  x-ray  peaks  found  on  the  Si-rich  side  of  Pd2Si  indicating 
both  Pd2Si  and  PdSi  to  be  present.  Later,  more  detailed  analysis  failed  to 
detect  any  PdSi  peaks  until  ~  50  At  %  SI.  Perhaps  unusually  high 

mobilities  of  Si  in  Pd2S1  near  room  temperature  led  to  the  indicated 


structural  rearrangement.  X-ray  analysis  of  the  deposited  films  at  the 
composition  of  Pd2Si  using  a  Read  camera  indicated  strong  texturing.  The 
particular  texture  is  that  expected  for  Pd2Si  growth  on  (111)  Si(Pd2Si 
surface  normal  is  approximately  9°  from  the  Si  surface  normal). 

For  the  elevated  temperature  (400°C)  phase  spread  deposition,  where 
the  entire  substrate  was  located  in  the  chamber  at  the  position  where  Pd2Si 
should  form,  the  entire  deposited  film  had  the  Pd2Si  structure  although 
electron  microprobe  analysis  showed  ~  6%  variation.  The  larger  region  of 
reaonsably  homogeneous  composition,  as  compared  to  the  lower  temperature 
phase  spread  deposition,  was  probably  due  to  the  increased  surface  mobility 
of  Pd  and  Si. 


4.  Electrical  Diagnostics  of  SiOw  Films 

The  designated  goal  of  this  research  is  to  determine  the  feasibility 

of  extracting  useful  and  reliable  data  from  non-contact  in-situ  monitoring 

of  Si 02  growth  by  the  A.C.  Impedance  technique.  The  proposed  experimental 

set-up  for  such  monitoring  Is  as  given  In  Fig.  6.  Ideally,  one  would 

oxidize  the  Si  samples  at  temperature  T  and,  at  periodic  time  intervals, 

one  would  perform  an  automated  measurement  from  a  desk  top  computer.  An 

Individual  measurement  would  consist  of  various  Impedance  readings 

(B,G,X,R,e '  ,e",C  R  ,  etc.).  When  optimized,  this  measurement  cycle 
P  P 

should  take  ~  1  mln-30  min  depending  on  the  need  for  lower  frequencies  (1 
mHz  -  1  Hz).  The  result  of  this  measurement  cycle  would  be  data  on  a 

floppy  disc  or  magnetic  tape  to  be  eventually  plotted  as  X-R,  B-G  and 

e'-e"  plots.  These  plots  would  contain  Information  on  the  combined 
Impedance  of  the  entire  sample  geometry.  One  would  hope  to  Isolate  and 
extract  the  resistance,  Rnv,  and  the  capacitance,  Cnw,  of  the  bulk  SiO?  and 


relate  this  to  the  thickness  of  the  growing  oxide,  tox,  via 

R  =  p(T)t  /A  and/or  C  =  e  e  A/t  .  This  process,  with  some  on-line 
ox  '  '  ox  ox  o  ox 

analysis  of  the  data  by  the  desk  top  computer,  should  yield  real-time 
monitoring  of  the  thickness  without  the  need  for  plotting  the  non-impedance 
data. 

The  success  of  this  experiment  depends  on  the  ability  to  view  the 

plotted  data  (X-R,  B-G,  e’-e")  and  sensitively  distinguish  one  component 

of  the  electrical  system  from  another.  To  perceive  the  problem,  consider 

the  sample  electrical  equivalent  circuit  given  in  Fig.  7.  The  dominant 

Impedances  will  be  the  Si  02/Si  interface,  the  bulk  Si  02  and  the  air  gap 

leading  to  the  reduced  equivalent  circuit  of  Fig.  8. 

Some  preliminary  C-V,  I-V  and  AC  measurements  at  room  temperature 

showed  that  *  Cox  at  25°C.  Using  conductivity  data  from  Srivastava  et 

al^13),  we  find  that  P  *  5  x  1015  and  p  *  1010  at  25°C  abd  850°C 

respectively  so  that  R  *■  104  -  107  w  at  850°C  for  t  *  1000  A  and 

ox  ox 


A  =  10  cm2  -  10"2cm2. 


Assuming  that  e  is  independent  of  T, 


Cox  *  3.5  x  10”7F  to  3.5  x  10”*°F  (as  A  changes  in  this  range). 

The  air  gap  brings  in  many  complicated  considerations.  It  must  be 
wide  enough  to  allow  oxygen  to  diffuse  freely  to  the  oxide  surface  and  yet 
it  must  be  narrow  enough  to  give  a  small  value  of  Cag  and  be  experimentally 

controllable.  If  A  =  10  cm2,  d  =  lum  ♦  C  _  “  10"10F  but  if  d  =  10“6  cm 

ag 
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♦  C  “  10  F.  Since  we  now  have  some  idea  of  the  parameter  magnitudes  in 
ag 

Fig.  8,  we  can  now  synthesize  the  sample  plots  X-R,  B-G  and  c"-e‘  if  we 
neglect  C^.  The  X-R  plots  Illustrate  the  salient  features  present  in  all 
the  plots  so  we  shall  consider  only  these.  In  Figs.  9(a)  and  9(b),  we  note 
the  great  difference  In  the  X  versus  R  plots  as  Cag  is  increased  from  10"10 
F  to  10”®  F.  In  the  former  case,  it  is  not  clear  that  we  can  readily 
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separate  the  effects  of  the  two  components.  From  this  shape  of  curve  we 
can  conclude  that  Cag  is  of  the  same  order  as  Cox  otherwise  such  a  strong 
coupling  would  not  have  occurred.  In  Fig.  9(b),  the  two  segments  of  the 
circuit  are  uncoupled  and  we  can  readily  distinguish  Cag  from  CQX. 

From  the  foregoing,  we  can  conclude  that  the  effect  of  an  arbitrary 
interface  impedance  will  be  generally  much  easier  to  pick  out  of  the  data 
base  when  the  two  initial  impedance  effects  are  "displayed"  at  sufficiently 
different  frequency  ranges.  This  is  not  true  when  Cag  =  10"*°  F  which  is 
the  experimentally  more  realizable  value.  Thus,  we  must  conclude  that  the 
use  of  this  technique  for  sensitively  studying  unknown  systems  is  only  on 
the  edge  of  experimental  viability  for  the  general  case.  For  specific 
situations  when  one  can  arrange  Cflg  to  be  very  different  from  CQX,  the 
technique  will  have  great  utility  for  studying  interface  effects. 

B.  Theoretical  Results 

1.  Film  Formation  Modeling  via  Computer  Simulation 
a.  Parametric  Investigations 

Details  of  film  formation  require  a  knowledge  of  the  map  of 

stored  strain  energy  in  the  substrate-film  combination.  However,  the 

available  techniques  for  evaluating  the  stress  tensor  components  and  the 

elastic  constants  of  a  two  phase  system  In  the  Immediate  vicinity  of  the 

Interface  region  are  quite  Inadequate  so  we  decided  to  develop  a 

generalized  statistical  thermodynamic  approach  to  the  problem.  The  first 

result  Is  that  individual  stress  and  elastic  formula  have  been  developed 

for  the  Interface  layers  by  calculating  the  "local"  pressure  variations 

across  the  Interface.  The  basic  results  for  the  stress,  c>  ,  and  elastic 

ij 
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constants  C 


are 


1  3A 


1  3‘A 


V  3n..3n 


where 


ij  'k i 


A  =  -  kT  inQ 


Q  =  / “•  J dxNdpNexp[-H(x,N  PN)/kT] 


Here,  A  is  the  Helmholtz  free  energy,  n.^  is  the  Lagrangian  strain 
parameter  and  H  is  the  Hamiltonian  for  the  system.  Equations  (9)  and  (10) 
are  always  evaluated  at  zero  strain  but  more  general  forms  have  been 
derived  in  the  vicinity  of  zero  strain.  Explicit  analytical  forms  have 
been  calculated  for  the  strain  derivatives  using  both  two-body  (Lennard- 
Jones)  and  three-body  (Axelrod-Teller)  potential  energy  functions.  From 
this  work  we  find  that  (1)  the  entire  variation  from  one  stress  component 
or  elastic  constant  to  another  is  determined  by  the  general  symmetry  of  the 
two  and  three  particle  distribution  functions  so  that  the  result  can  be 
separated  Into  a  geometrical  combination  factor  and  a  material  combination 
factor,  (11)  new  and  unique  symmetry  patterns  are  introduced  via  the  three- 
body  potential  so  that,  In  principal,  all  21  elastic  constants  can  be 
unique  (rather  than  just  15  for  the  purely  two-body  PEF  case)  and  (111)  our 
equations  are  perfectly  general  and  apply  to  both  bulk  solid  as  well  as 
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non-homogeneous  regions  like  Interfaces. 

As  mentioned  In  the  previous  report(^),  the  adlayer-substrate 
interaction  effects  are  being  investigated  using  a  system  of  1000  particles 
(~  400  adlayer,  X,  and  ~  600  substrate,  Ag(lll))  using  only  two-body  forces 
for  a  silver  substrate  and  adlayers  of  Cu(lll),  Au(lll),  Ni ( 100  and 
Ag(lll).  The  bottom  three  layers  of  the  six  substrate  layers  are  held 
fixed  and  all  thermodyanic  quantities  were  calculated  considering  only  the 
“core"  region  of  the  hexagonal  slab  (to  eliminate  edge  effects).  The 
calculations  generated  potential  energy  profiles,  density  profiles,  radial 
distribution  function,  virial  profiles,  vibrational  spectra  and  stress 
components.  Because  of  its  differential  character  the  virial  reflects  the 
mismatch  of  the  sites  at  the  interface^);  however,  some  of  the  stress 
components  give  the  most  interesting  resuls. 

In  spite  of  the  11.5%  mismatch  in  the  Cu/Ag  system,  the  interface  is 

coherent  with  increased  compressive  forces  in  the  substrate  parallel  to  the 

interface  while  the  adlayer  itself  is  in  the  tensile  state.  For 

°xx  *n  the  only  the  top  Cu  layer  shows  a  stress  dipole 

behavior  while,  for  a ^  (perpendicular  to  the  Interface)  a  much  larger 

dipole  effect  exists  for  the  first  layer  with  a  smaller  dipole  effect 

existing  for  the  lower  layers  and,  with  o  (In-plane  shear  stress),  no 

xy 

stress  dipole  effect  exists  at  all.  The  results  for  o ^  and  are 

given  in  Fig.  10. 

For  the  Au/Ag  case,  although  the  atomic  mismatch  is  only  0.1%,  a  very 
strong  stress  effect  is  observed  because  the  Au-Au  interaction  energy  is 

much  stronger  than  the  Au-Ag  interaction  energy.  The  vibrational  frequency 
spectra  for  the  Ag  atoms  at  the  Au/Ag  and  Cu/Ag  interfaces  also  displayed 
considerable  differences  in  peak  positions  and  shapes.  The  largest  effect 


was  found  for  the  z-component  of  the  vibrational  spectra  due  to  the  heavy 
mass  of  the  Au  atoms. 

For  the  Ni/Ag  case,  the  interface  exhibits  a  structural  channel 

direction  with  a  very  different  stress  picture  along  versus  normal  to  the 

channel  direction.  Along  the  channel  direction  there  is  no  strong  coupling 

of  atomic  positions  across  the  interface  so  there  is  negligible  stress 

coupling,  however,  perpendicular  to  the  channel  direction,  strong 

correlation  of  atomic  positions  and  strong  stress  coupling  exists.  This  is 

demonstrated  in  the  comparison  of  o  and  a  presented  in  Fig.  11. 

xx  yy 

We  find  that  our  method  of  resolving  the  stress  contributions  by 
atomic  layers  is  a  workable  technique  with  real  merit  because  the 
calculated  results  can  be  understood  in  terms  of  atomic  processes.  It 
looks  as  If  (i)  the  stress  dipole  effect  may  be  useful  as  a  yardstick  for 
measuring  the  surface  influence  penetration  in*.o  the  bulk,  (ii)  the  cross¬ 
interaction  parameters  (ex.^g)  an  important  guide  to  anticipate 

tangential  shear  strain  dependence  and  (iii)  stress  effects,  in  the  absence 
of  atomic  mismatch,  may  provide  information  on  clustering  tendencies  and 
process  activation  barriers.  Extension  of  this  work  to  include  three-body 
forces  and  periodic  boundary  conditions  seems  to  be  an  important  step  to 
take. 


(b)  Simulation  of  Systems  with  Si  and  0  Atoms 


A  Monte-Carlo  procedure,  based  on  two-  and  three-body  potentials 
was  employed  to  simulate  crystalline  o-cristobalite  and  a-quartz 
structures  at  T  =  50°K  and  T  *  300°K.  Iterations  were  carried  out  beyond 
the  complete  equilibration  stage  (monitored  by  variation  of  the  total 
energy).  Both  the  radial  distribution  function  (RDF)  and  the  angular 


distribution  function  (ADF)  for  Si -0-Si  and  0-Si -0  are  determined  and  are 
shown  in  Figs.  12  and  13  for  o-cristobalite  at  300°K.  Analysis  of  these 
results  Indicates  that  the  initial  crystalline  structures  for  both  phases 
are  quite  stable  implying  that  the  potential  energy  parameters  used  in 
these  computations  provide  deep  and  stable  minima  at  the  temperatures 
considered.  No  changes  are  found  in  the  RDF  and  ADF  when  the  number  of 
atoms  was  increased  from  96  to  216.  The  data  of  Fig.  14  for  Si-O-Si  angle 
peaks  at  ~  155°  which  is  a  little  high  compared  to  the  experimental  data. 

At  T  >  2000°K,  both  the  a-quartz  and  the  a-cristobalite  displayed 
liquid-like  characteristics  in  the  RDF's  as  expected.  This  is  an  important 
result  because  the  Born-Mayer-Huggins  PEF  used  by  almost  all  other 
investigators  for  Si O2  requires  them  to  go  to  T  ~  5000°K  before  liquid- 
like  characteristics  appear. 

Two  isolated  clusters,  Si O4  and  Si207,  were  treated  separately  with 
the  Monte-Carlo  procedure  at  T  *  300°K.  Iterations  of  long  duration 
indicated  that  clusters  initially  generated  as  tetrahedral  units  preserved 
their  tetrahedral  structure. 

We  next  attempted  to  simulate  the  condition  of  the  Si 02  surface  using 
the  Monte-Carlo  technique.  Figure  14(a)  shows  an  unrelaxed  (lOO)surface  of 
o-cristobalite  while  Fig.  14(b)  shows  the  reconstructed  surface 
condition.  In  this  computation,  the  top  four  layers  of  the  system  were 
allowed  to  relax  while  the  remainder  of  the  lattice  was  held  fixed. 
Periodic  boundary  conditions  were  used  on  our  basic  system  of  72  Si  and  144 

0  atoms  (3x3x3x12  atoms/unit  cell).  Comparing  Figs.  14(a)  and  (b),  we  note 

that  oxygen  atoms  of  the  first  layer  tend  to  form  pairs  with  some 

regularity  while  0  atoms  of  the  third  layer  do  not  move.  These  pairs 
appear  to  be  in  some  sense  correlated.  We  have  not  confirmed  any  long 
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*a 


Number  Substrate 


111  Si 


111  Si 


TABLE  I 

SiC  Sputter  Depositions 
Clean*  T,  SiC,  Sn 


Purpose 


as  rec'd 

900 °C,  600W,  $ 

heater  cycle 

Hf,  DI 

II  ll 

dep 

as  rec'd 

25° C,  4>  »  unstable 

Sn  rate  test 

MeOH,  DI 

25°  C,  <J>  ,  180W 

ii 

li 

108W 

li 

li 

44W 

ll 

MeOH,  DI 

25°  C,  600W,  4> 

SiC  rate  test 

A,HF,MeOH,DI 

900° C,  <J>  ,  4> 

heater  cycle 

II 

900“ C,  600W,  166W 

dep  with  Sn 

A,  DI 

"  "  156W 

"  on  sapphire 

C,  DI 

*©- 

heater  c>>ctr. 

500 °C,  600W,  110W  new  T 
600 “C,  600W,  108W 
700°C,  600W,  101W 
800°C,  600W,  104W 
"  "  96W  new  sn 


TABLE  II 


SiC  Sputter  Depositions 

A.  600  W,  2"  system,  6  l/2"source-substrate. 


Number 

Substrate 

T 

Tin  Power 

Feature  Size  (ym) 

82-617 

(lll)Si 

900"  C 

Heater  Cycle 

3. 5-4.0 

82-619 

II 

750°C 

- 

0.6-1. 2 

82-620 

H 

750°C 

20  W 

0.6-1. 2,  1.3-2. 3 

82-621 

n 

750°C 

50  W 

1.0-2. 6 

82-622 

ii 

900  #C 

20  W 

1. 0-2.0,  1.3-3. 8,  4. 2-6. 5 

82-623 

n 

900  “C 

50  W 

0.6-1. 3,  1.3-3. 2,  3. 2-4. 5 

B.  500 

Number 

W,  3"  system,  3" 

Substrate 

source- 

T 

substrate 

Relative  X-ray 

Peak  Height 

Reflections  in  Read 

Pattern 

82-123 

0x(100)Si 

900°C 

0.4 

111,  220,  222 

84-124 

(oooi)ai2o3 

800°C 

100.0 

111 

82-125 

MgO  A1203 

800°  C 

3.6 

111,  200,  220,  222,  1010 

81-233* 

(lll)Si 

750°C 

11.5 

111 

★ 

Best  previous  deposition 


TABLE  III 


t£< 


* 


Pd-Si  Phase  Spread  Deposition 

Si  power  =  570  W  RF  -*■  272  A/min  ;  P  =  3  ym 

ar 

Pd  current  =  0.11  A  -*•  200  A/min  ;  t  =  93  min 


Composition 

Structure 

At  %  Si 

At  %  Pd 

0 

100 

2 

98 

4 

96 

Pd 

6 

94 

8 

92 

10 

90 

14 

86 

Amorphous  &  Pd 

17 

83 

Amorphous 

19 

81 

II 

23 

77 

Amorphous  &  Pd^Si 

33 

67 

Pd2Si 

36 

64 

41 

59 

Pd2Si  &  PdSi 

44 

56 

52 

48 

PdSi 

66 

34 

? 

76 

24 

? 
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FIGURE  CAPTIONS 


1.  Oxygen  concentration  for  reactively  synthesized  films  is  shown  as  a 

function  of  the  ratio  of  the  oxygen  (O2 )  and  silicon  incidence  rates 
for  Type  IIA  and  I  IB  experiments  and  for  the  thermal  source  results  of 
Ritter^6). 

2.  Oxygen  concentration  (x  in  Si 0X )  is  shown  as  a  function  of  dynamic 

exposure,  Lp,  for  films  deposited  using  the  Type  11A  and  I IB 
experimental  arrangements  and  the  thermal  data  of  Ritter^. 

3.  Experimental  results  for  Type  IIA  (I)  and  Type  I  IB  (#,H,A)  are 

compared  to  the  reported  results  of  oxygen  coverage  on  Si (111) 
surfaces  as  determined  by  HRES  (□)  ellipsometry  (O)  and  Auger 
analysis  (A)  assuming  that  dynamic  exposure  (Lq)  and  exposure  (L)  are 
equivalent. 

4.  X  in  Si0x  is  shown  as  a  linear  function  of  dynamic  exposure,  Lp, 

(Torr-sec)  allowing  estimation  of  the  sticking  coefficients,  S,  for 
oxygen  in  the  Type  IIA  and  Type  I IB  experiments  (see  Eq.  (1)). 

5.  The  intramolecular  potential,  Vp,  is  plotted  versus  distance 
illustrating  the  potential  wells  for  both  chemisorption  and 
physisorption  Is  schematically  shown.  The  physisorbed  state  (DEF) 
represents  the  mobile  precursor  with  an  adsorption  energy  qp.  If  the 
activation  energy,  Ep,  for  transition  from  the  precursor  to  the 
chemisorbed  state  (ABC)  Is  less  than  qp,  the  activation  energy  for 
chemisorption  Is  zero. 

6.  Schematic  Illustration  of  experimental  set-up  for  non-contact,  in-situ 
monitoring. 

7.  Sample  electrical  equivalent  circuit  for  Fig.  6. 


8.  Reduced  equivalent  circuit  for  Fig.  6. 

9.  Plots  of  circuit  reactance,  x,  versus  circuit  resistance,  R,  with  Cox 
=  5  x  10"11  F  and  RQX  =  2  x  105  ft  at  a  range  of  frequencies  for  the 


two  cases  (a)  Cag 

=  10‘10F  and  (b)  Cag 

U. 

00 

o 

M 

• 

10.  Stress  components 

0 

XX 

and  ozz  for 

the  Cu 

(111) 

adlayer/Ag 

(111) 

substrate  system 

as  a 

function  of 

atomi c 

position  normal  to  the 

interface. 

11.  Stress  components 

°xx 

and  ozz  for 

the  Ni 

(100) 

adlayer/Ag 

dll) 

substrate  system 

as  a 

function  of 

atomi c 

position  normal  to  the 

Interface. 

12.  Radial  Distribution  Function  for 

a-cristobalite 

(300°K) 

as  a 

function  of  bond 

length:  (a)  Si-Si 

bonds. 

(b)  0-0 

bonds,  (c) 

Si-0 

bonds,  (d)  total  bonds. 

13.  Angular  Distribution  Function  for  a-cristobalite  (300°K)  as  a 
function  of  angle  for  both  Si -0-Si  and  O-Si-O. 

14.  (100)  Surface  of  a-cristobalite 

(a)  before  relaxation 

(b)  after  relaxation. 
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